Search for: All records

Abstract Specialist insects are especially susceptible to loss of genetic diversity in the face of habitat destruction and fragmentation. Implementing effective conservation practices for specialist insects will benefit from knowledge of population structure and genetic diversity. Because insects are hyper-diverse, characterizing the population structure of all species within the insect community is untenable, even if focused within a particular habitat type. Thus, concentrating on a single species specialized to a particular habitat type is needed to infer general trends. Here, we investigate the range-wide population genetics ofTetraopes texanusHorn 1878 (Coleoptera: Cerambycidae), which provides a useful model of grassland insects due to its’ habitat specificity and unique biology.Tetraopes texanusoccurs primarily in Texas and Oklahoma, into Northern Mexico, and possibly into eastern New Mexico but also occurs in Black Belt prairies of Mississippi and Alabama. Mitochondrial and nuclear DNA (RAD-seq) analysis identified two distinct population clusters ofT. texanuscorresponding to the Texas and Oklahoma population and the Mississippi and Alabama population. Demographic models indicate ongoing, though incomplete, isolation of the two populations, with estimated dates of divergence in the mid-Pleistocene, coinciding with the end of a glacial period and a shift in glacial interval. These results can inform conservation of grassland adapted insects and offers insight to the biogeography of the Gulf Coastal Plain. 

Abstract Bumble bees (Bombus) exhibit exceptional diversity in setal body color patterns, largely as a result of convergence onto multiple Mullerian mimicry patterns globally. When multiple species cross the same sets of mimicry complexes, they can acquire the same color polymorphisms, providing replicates of phenotypic evolution. This study examines the genetic basis of parallel color pattern acquisition in three bumble bee taxon pairs in western North America that shift between orange-red and black mid-abdominal segmental coloration in Rocky Mountain and Pacific Coastal mimicry regions: polymorphic Bombus vancouverensis and B. melanopygus, and sister species B. huntii and B. vosnesenskii. Initial gene targets are identified using a genome-wide association study, while cross-developmental transcriptomics reveals genetic pathways leading to final pigmentation genes. The data show all three lineages independently target the regulatory region of a segmental-fate determining Hox gene, Abdominal B (Abd-B), for this color transition. For B. vancouverensis and B. melanopygus, this involves different deletions in the same location, and all mimicry pairs differentially express Abd-B and ncRNAs in this locus. Transcriptomics reveals a shared core gene network across species, where Abd-B interacts with nubbin and pigment enzyme ebony to decrease black melanin production in favor of paler, redder morphs. Expression of multiple genes in the melanin biosynthesis pathway is modified to promote this phenotype, with differing roles by taxon. Replicated morphologies unveil key genes and a Hox gene hotspot, while enabling evolutionary tracking of genetic changes to phenotypic changes and informing how gene regulatory networks evolve. 

Abstract We present the first chromosome-level genome assembly for Bombus pensylvanicus, a historically widespread native pollinator species that was distributed across eastern North America but has subsequently undergone declines in range area and local relative abundance. This species has been of significant interest as a model for understanding both patterns and possible causes of bumble bee decline in the region, including the role of genetic variation. Here we present a chromosome-level reference genome assembled using Pacific Biosciences singe-molecule HiFi sequences and Hi-C data and annotated using evidence derived from RNA sequencing of multiple tissue types. The B. pensylvanicus genome has a total length of ∼352.6 Mb and was assembled into a total of 224 scaffolds, with 19 primary pseudomolecules representing putative chromosomes and an N50 = 14.872 Mb. Annotation with the Eukaryotic Genome Annotation Pipeline—External (EGAPx) identified 11,411 genes (10,263 protein coding), and BUSCO analysis of 5,991 Hymenoptera-specific BUSCO groups indicated a completeness for the proteins of 99.0% (98.6% single-copy, 0.5% duplicated) and for the genome of 98.5% (98.2% single-copy, 0.3% duplicated). We present synteny analyses with other recently assembled Bombus genomes representing different subgenera and examine the distribution of repetitive regions of the genome relative to the distribution of genes and noncoding RNAs. 

Abstract The Hunt bumble bee, Bombus huntii, is a widely distributed pollinator in western North America. The species produces large colony sizes in captive rearing conditions, experiences low parasite and pathogen loads, and has been demonstrated to be an effective pollinator of tomatoes grown in controlled environment agriculture systems. These desirable traits have galvanized producer efforts to develop commercial Bombus huntii colonies for growers to deliver pollination services to crops. To better understand Bombus huntii biology and support population genetic studies and breeding decisions, we sequenced and assembled the Bombus huntii genome from a single haploid male. High-fidelity sequencing of the entire genome using PacBio, along with HiC sequencing, led to a comprehensive contig assembly of high continuity. This assembly was further organized into a chromosomal arrangement, successfully identifying 18 chromosomes spread across the 317.4 Mb assembly with a BUSCO score indicating 97.6% completeness. Synteny analysis demonstrates shared chromosome number (n = 18) with Bombus terrestris, a species belonging to a different subgenus, matching the expectation that presence of 18 haploid chromosomes is an ancestral trait at least between the subgenera Pyrobombus and Bombus sensu stricto. In conclusion, the assembly outcome, alongside the minimal tissue sampled destructively, showcases efficient techniques for producing a comprehensive, highly contiguous genome. 

Abstract Understanding the myriad avenues through which spatial and environmental factors shape evolution is a major focus in biological research. From a molecular perspective, much work has been focused on genomic sequence variation; however, recently there has been increased interest in how epigenetic variation may be shaped by different variables across the landscape. DNA methylation has been of particular interest given that it is dynamic and can alter gene expression, potentially offering a path for a rapid response to environmental change. We utilized whole genome enzymatic methyl sequencing to evaluate the distribution of CpG methylation across the genome and to analyze patterns of spatial and environmental association in the methylomes of two broadly distributed montane bumble bees (Bombus vancouverensis Cresson and Bombus vosnesenskii Radoszkowski) across elevational gradients in the western US. Methylation patterns in both species are similar at the genomic scale with ∼1% of CpGs being methylated and most methylation being found in exons. At the landscape scale, neither species exhibited strong spatial or population structuring in patterns of methylation, although some weak relationships between methylation and distance or environmental variables were detected. Differential methylation analysis suggests a stronger environment association in B. vancouverensis given the larger number of differentially methylated CpG's compared to B. vosnesenskii. We also observed only a handful of genes with both differentially methylated CpGs and previously detected environmentally associated outlier SNPs. Overall results reveal a weak but present pattern in variation in methylation over the landscape in both species. 

Abstract Unraveling molecular mechanisms of adaptation to complex environments is crucial to understanding tolerance of abiotic pressures and responses to climatic change. Epigenetic variation is increasingly recognized as a mechanism that can facilitate rapid responses to changing environmental cues. To investigate variation in genetic and epigenetic diversity at spatial and thermal extremes, we use whole genome and methylome sequencing to generate a high-resolution map of DNA methylation in the bumble bee Bombus vosnesenskii . We sample two populations representing spatial and environmental range extremes (a warm southern low-elevation site and a cold northern high-elevation site) previously shown to exhibit differences in thermal tolerance and determine positions in the genome that are consistently and variably methylated across samples. Bisulfite sequencing reveals methylation characteristics similar to other arthropods, with low global CpG methylation but high methylation concentrated in gene bodies and in genome regions with low nucleotide diversity. Differentially methylated sites (n = 2066) were largely hypomethylated in the northern high-elevation population but not related to local sequence differentiation. The concentration of methylated and differentially methylated sites in exons and putative promoter regions suggests a possible role in gene regulation, and this high-resolution analysis of intraspecific epigenetic variation in wild Bombus suggests that the function of methylation in niche adaptation would be worth further investigation. 

Abstract The gut microbiome is influenced by host species and the environment, but how the environment influences the microbiome of animals introduced into a new ecosystem has rarely been investigated. Freshwater mussels are aquatic fauna, with some threatened or endangered species propagated in hatcheries and introduced into natural systems as part of conservation efforts. The effects of the environment on the freshwater mussel gut microbiome were assessed for two hatchery-propagated species (Lampsilis ovata, Lampsilis ornata) introduced into rivers within their natural range. Mussels were placed in rivers for 8 weeks, after which one subset was collected, another subset remained in that river, and a third subset was reciprocally transplanted to another river in the same river basin for a further 8 weeks. Gut microbiome composition and diversity were characterized for all mussels. After the initial 8 weeks, mussels showed increased gut bacterial species richness and distinct community composition compared to hatchery mussels, but gut microbiome diversity then decreased for mussels that remained in the same river for all 16 weeks. The gut bacterial community of mussels transplanted between rivers shifted to resemble that of mussels placed initially into the recipient river and that remained there for the whole study. All mussels showed high proportions of Firmicutes in their gut microbiome after 8 weeks, suggesting an essential role of this phylum in the gut of Lampsilis species. These findings show that the mussel gut microbiome shifts in response to new environments and provide insights into conservation strategies that involve species reintroductions. 

ABSTRACT Bumble bees are common in cooler climates and many species likely experience periodic exposure to very cold temperatures, but little is known about the temporal dynamics of cold response mechanisms following chill exposure, especially how persistent effects of cold exposure may facilitate tolerance of future events. To investigate molecular processes involved in the temporal response by bumble bees to acute cold exposure, we compared mRNA transcript abundance in Bombus impatiens workers exposed to 0°C for 75 min (inducing chill coma) and control bees maintained at a constant ambient temperature (28°C). We sequenced the 3′ end of mRNA transcripts (TagSeq) to quantify gene expression in thoracic tissue of bees at several time points (0, 10, 30, 120 and 720 min) following cold exposure. Significant differences from control bees were only detectable within 30 min after the treatment, with most occurring at the 10 min recovery time point. Genes associated with gluconeogenesis and glycolysis were most notably upregulated, while genes related to lipid and purine metabolism were downregulated. The observed patterns of expression indicate a rapid recovery after chill coma, suggesting an acute differential transcriptional response during recovery from chill coma and return to baseline expression levels within an hour, with no long-term gene expression markers of this cold exposure. Our work highlights the functions and pathways important for acute cold recovery, provides an estimated time frame for recovery from cold exposure in bumble bees, and suggests that cold hardening may be less important for these heterothermic insects. 

Abstract Studies of species that experience environmental heterogeneity across their distributions have become an important tool for understanding mechanisms of adaptation and predicting responses to climate change. We examine population structure, demographic history and environmentally associated genomic variation in Bombus vosnesenskii , a common bumble bee in the western USA, using whole genome resequencing of populations distributed across a broad range of latitudes and elevations. We find that B. vosnesenskii exhibits minimal population structure and weak isolation by distance, confirming results from previous studies using other molecular marker types. Similarly, demographic analyses with Sequentially Markovian Coalescent models suggest that minimal population structure may have persisted since the last interglacial period, with genomes from different parts of the species range showing similar historical effective population size trajectories and relatively small fluctuations through time. Redundancy analysis revealed a small amount of genomic variation explained by bioclimatic variables. Environmental association analysis with latent factor mixed modelling (LFMM2) identified few outlier loci that were sparsely distributed throughout the genome and although a few putative signatures of selective sweeps were identified, none encompassed particularly large numbers of loci. Some outlier loci were in genes with known regulatory relationships, suggesting the possibility of weak selection, although compared with other species examined with similar approaches, evidence for extensive local adaptation signatures in the genome was relatively weak. Overall, results indicate B. vosnesenskii is an example of a generalist with a high degree of flexibility in its environmental requirements that may ultimately benefit the species under periods of climate change.